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We investigate the violation of local realism in Bell tests involving homodyne measurements 
performed on multimode continuous-variable states. By binning the measurement outcomes in an 
appropriate way, we prove that the Mermin-Klyshko inequality can be violated by an amount that 
grows exponentially with the number of modes. Furthermore, the maximum violation allowed by 
quantum mechanics can be attained for any number of modes, albeit requiring a quantum state 
that is rather unrealistic. Interestingly, this exponential increase of the violation holds true even for 
simpler states, such as multipartite GHZ states. The resulting benefit of using more modes is shown 
to be significant in practical multipartite Bell tests by analyzing the increase of the robustness 
to noise with the number of modes. In view of the high efficiency achievable with homodyne 
detection, our results thus open a possible way to feasible loophole-free Bell tests that are robust to 
experimental imperfections. We provide an explicit example of a three-mode state (a superposition 
of coherent states) which results in a significantly high violation of the Mermin-Klyshko inequality 
(around 10%) with homodyne measurements. 



I. INTRODUCTION 

The non-local nature of quantum mechanics is unar- 
guably one of its most counterintuitive aspects, mark- 
ing a fundamental departure from the classical picture of 
physical systems. First recognized by Einstein, Podolsky, 
and Rosen [l[ , the intrinsic non- locality of the quantum 
theory was put into a testable form via the celebrated 
Bell inequalities 0). Since the first experimental viola- 
tion of a Bell inequality *3| i a great variety of tests have 
been performed. In all of them, a local realistic model 
is shown to be incompatible with the actual experimen- 
tal observations. Nonetheless, a fully conclusive evidence 
for rejecting local realism has not been achieved yet, for 
all the experiments performed so far suffer from loop- 
holes, forcing us to rely on supplementary assumptions 
in order to reject local realism. Specifically, a conclusive 
experiment should close two main loopholes: the locality 
loophole (namely, the measured correlations must be col- 
lected with spacelike separated events) and the detection- 
efficiency loophole (namely, the proportion between de- 
tected and undetected events has to be high enough for 
the data to be fully representative of the whole ensemble 
without the need to assume a "fair sampling"). 

There is a large variety of quantum systems for which 
a test of local realism may be envisaged. However, the 
quest for a loophole- free Bell test has recently focused the 
research towards experiments involving propagatinglight 
modes measured with homodyne detectors [J, SB 13 • 
The advantage of such a "continuous-variable" approach 
is twofold. First, light modes can be easily sent to space- 
like separated detectors suffering only a tolerable degree 
of decoherence. Second, the current technology of homo- 
dyne detectors achieves a degree of detection efficiency 
high enough to potentially close the detection efficiency 
loophole. On the other hand, such an approach involves 
also drawbacks whose resolution is still challenging. The 
main issue comes from the fact that the results of ho- 



modyne measurements can be described by means of the 
Wigner function. Thus, in order to avoid a local hid- 
den variable description of the measured correlations, 
one should necessarily perform the test with a state en- 
dowed with a non-positive Wigner function (otherwise, 
the Wigner function is a genuine probability distribu- 
tion, which provides an explicit local realistic model of 
the data). Even if single-mode states of traveling light 
have already been generated with a non-positive Wigner 
function [1, §, such (non-Gaussian) states are hard to 
master in the laboratory. Among the above mentioned 
proposals of Bell tests with homodyne detection, those 
nearer to an experimental realization only give a small 
violation of Bell inequalities [y, [7[ , whereas higher viola- 
tions involve states whose actual generation seems rather 
unrealistic [qj. In this context, the search for Bell tests 
involving feasible resources and giving, at the same time, 
violations that are high enough to be robust against ex- 
perimental noise is very desirable. This is the thrust of 
the present paper. 

Bell tests relying on homodyne detection typically 
involve the following scenario. Two parties perform 
spacelike separated homodyne measurements by ran- 
domly choosing each between two settings, thus mea- 
suring two quadratures of the incoming electromagnetic 
fields. The violation of the Bell Clauser-Horne-Shimony- 
Holt (CHSH) inequality ^ is then tested. Since this 
inequality is devised for two-outcome measurements, the 
collected data, which are distributed (ideally) in a con- 
tinuous way, have to be discretized; such a procedure is 
referred to as binning process. In this work, we will gen- 
eralize this scenario to more than two parties, consider- 
ing the Mermin-Klyshko (MK) type of Bell inequalities, 
which involves m p arties, two measurement settings, and 
two outcomes [11] . Our main motivation is that higher 
violations of local realism are expected as the number 
of parties increases [l2j], which may help in the search 
for feasible proposals of loophole-free Bell tests able to 
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tolerate the experimental noise. 

The earlier investigations of Bell tests in infinite- 
dimensional spaces with multipartite settings have con- 
sidered several scenarios. In Ref. [IH, a test based on 
the measurement of the light field parity was envisaged. 
It was found that the violation does not increase expo- 
nentially, as one would have expected, but this may be 
due to the fact that no optimization over the possible 
measurement settings was performed [l4| . In Ref. [l5j in- 
stead, a maximal violation of Mermin-Klyshko inequality 
was found for continuous variable states, considering the 
measurement of a different class of operators which can 
be seen as the continuous-variable analogue of the spin 
operators. However, both these approaches deal with 
non-Gaussian measurements described by a non-positive 
Wigner function, which are far from the reach of cur- 
rent detection technology. Recently, another approach 
has been introduced that does not rely on the use of in- 
equalities with discrete outcomes, thus avoiding the need 
of a binning procedure [161 ] . There, quadrature measure- 
ments via homodyne detection are considered, and the 
measurement outcomes are used to directly test a novel 
Bell inequality for continuous variables. Interestingly, the 
authors find that the violation of the local realistic bound 
is exponential in the number of parties. However, a possi- 
ble experimental implementation is still very challenging 
as it would require at least ten space-like separated ho- 
modyne measurements. It is then unclear whether such a 
novel approach can give advantages, from a practical per- 
spective, over non-locality tests involving binning strate- 
gies. 

In this paper, we will stick with the use of discrete- 
variable Bell inequalities in the tests of non-locality for 
continuous- variable states, thus using some binning pro- 
cedure. We will address the following question: Is it 
possible to have an exponential increase of the violation 
of local realism in a test involving quadrature measure- 
ments of m modes and considering the Mermin-Klyshko 
inequalities? We will answer this question by the af- 
firmative, providing specific examples of states exhibit- 
ing such a behavior. These states belong to the class 
of photon- number correlated states, and generalize to an 
arbitrary number of modes the approach of Ref. [4J . Re- 
markably, we will show that it is even possible to reach 
the maximal violation of the Mermin-Klyshko inequali- 
ties for to modes by properly choosing the binning pro- 
cedure and the quantum state, thereby extending the re- 
sults obtained in Ref. [f| for two modes. 

The theoretical interest of our results resides in the 
fact that we are not considering a direct mapping from 
continuous- to discrete-variable states. Thus, the possi- 
bility to obtain a maximal violation of Mermin-Klyshko 
inequality for any number of parties is a priori non- 
trivial. Furthermore, we may have anticipated that the 
binning procedure causes an irreparable loss of informa- 
tion, preventing the possibility to reach a maximum vi- 
olation. From an experimental perspective, our results 
imply that various strategies allow for an increase of the 



violation of locality in homodyne-based Bell tests. This, 
in turn, gives specific insights in the search of multipar- 
tite states appropriate for an actual experiment. 

The paper is organized as follows. In Sec. [Til we ana- 
lyze a Bell test involving multimode states with perfect 
correlations in the number of photons in all modes. We 
analytically show that the violation of local realism can 
grow exponentially with the number of modes involved. 
This is true despite the fact that a simple sign binning 
strategy is used. Nonetheless, we see that the resulting 
violation does not reach the maximum value achievable 
within quantum mechanics. However, we show in Sec. IIIII 
that a different binning strategy, properly tailored for 
some specific states, allows attaining the maximal vio- 
lation of locality achievable within quantum mechanics. 
Then, in Sec. HVl we address the robustness of these tests 
of nonlocality. We see that, at least for the analyzed noise 
model, the tolerable noise increases with the number of 
modes, making the violation better testable in practice. 
In Sec. El we exploit all these results from a practical 
oriented perspective. In particular, we give an explicit 
example of a class of three-mode states which, on the 
one hand, exhibit a gain in the violation by considering 
more than two modes, and, on the other hand, may be 
implemented with the near future technology. Finally, we 
close the paper in Sec lVIl with some concluding remarks. 

II. CORRELATED PHOTON NUMBER STATES 

We show in this Section that in order to obtain an ex- 
ponential violation of local realism with quadrature mea- 
surements, it is possible to rely on a simple binning strat- 
egy and the standard Mermin-Klyshko Bell inequalities. 
In achieving this goal we generalize to m parties the work 
of Ref. [3] , where the bipartite scenario was analyzed. 

Before proceeding, let us recall the general form of 
the Mermin-Klyshko Bell inequalities. Consider two di- 
chotomic observables Ot and 0' t for each party t. The 
Mermin-Klyshko inequalities are based on the recursive 
definition of the Bell operator 

B t = ^ ® (O t + 0[) + ^ ® (O t - 0' t ), (1) 

where Bx = 2 0i, B[ = 2 0[, and B' t denotes the same 
expression as B t but with all Ot's and O^'s exchanged 
[12|. The Mermin-Klyshko inequality for m parties then 
reads 

B m = \(B m }\<2. (2) 

In the case of two parties (m = 2), this expression reduces 
to the well-known CHSH inequality, for which 

B 2 = O x % 2 + O x % 0' 2 + 0[ ® 2 - 0[ ® 0' 2 (3) 

Let us now consider a generic photon-number corre- 
lated state of m bosonic modes 

oo 

l*> = X>|r) 1 |r> 2 ...|r) tnl (4) 

r=0 
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With J2r \°r 



1. The two observables Ot and Oi to one obtains that 



be measured are chosen as the two quadratures of the 
electromagnetic field X(9t) and X(6' t ) (corresponding to 
angles 8 and 8'). These can be measured by applying an 
homodyne detection on each mode t. In the following, 
we use the notation Xt for the quadrature X{6t), and 
x t for the outcome of its measurement (and similarly for 
primed quantities). 

The joint probability to obtain the outcomes x%, x m 
when measuring the quadratures X\, ...,X m is given by 

oo 

V(x- L ,...,X rn ) = 1 1 <rCl | ... m {x m \ *) | 2 = ^ ° r °* 

r,s=0 

P i<t>(r-s) ™_ 

]Je- x tH r {x t )H s (x t ) , 
t=i 

(5) 



{-K2 r+S r\s\) m / 2 



where <fi = 0± + - ■ ■ 8 m , H t (x) is the Hermite polynomial of 
degree t, and \xt) are the eigenvectors of the quadrature 
operator X t . 

Consider the following simple binning strategy, known 
as sign binning: when the result of a quadrature measure- 
ment falls in the domain M.^ the value ±1 is associated 
to it. Now, let us start by calculating the probability 
V+1,.,.,+1 that a result is observed in all the mea- 

suring sites: 

/>oo y* oo 

V+i +i = / dxi--- dx m V(xi,...,x m ) 

Jo Jo 

oo 

= E c 



i(j>(r-s) 



r,s—0 



n 



_-, JO 



t=l 



(7r2''+ s r!s!) m / 2 
dx t e- x tH r (xt)H s (xt). (6) 



The integrals above can be evaluated by recalling the 
following properties of Hermite polynomials for r =/= s, 

dxe- x2 H r (x)H s {x) = ^—l[F{r, S )-F{s,r)}, 

r-s 

(7) 



where we defined F(r, s) as 



(8) 



with r being the Gamma function. For r — s, one has 
instead 



poo 

/ dxe- x2 H*(x)=2 r - 1 r\V^, 
Jo 

Defining the functions 

0(0, m) = 2^2Re(c r c*) g r , s ((f>,m) , 



(9) 



,(<j>,m) 



7r2''+ s 



l/2 r 



r'.s 

x cos[cj)(r — s)] 



F(r,s) - F(s,r) 



(10) 



G{<j>,m) . 



(11) 



The other probabilities can be obtained in a similar way. 
Let us define the multi- index d = (d\, ■ ■ ■ d m ), with 
d t = ±1 denoting the measurement outcome obtained for 
mode t after binning. Then, the joint probability for a 
generic collection d of measurement outcomes will be in- 
dicated by Vd- It can be obtained similarly as V+i^.^+i 
by recalling that an Hermite polynomial of even (odd) 
degree is an even (odd) function, namely 



1 



a{d)Q{<t>,m) 



(12) 



where cr(d) = ]^L =1 d t . Now, we are in the position to cal- 
culate the generic correlation function between the mea- 
surement results E((j), m) (note that it only depends on 
the sum of the angles 4>). By definition, we have 



£(</>, m) =^a(d)7> d , 

d 



(13) 



where the sum goes over all possible collections d of mea- 
surement outcomes. Since the number of possible col- 
lections for which <r(d) = 1 is equal to that for which 
tr(d) = — 1, one finally has, by substituting Eq. (fT2| into 
Eq. that 



E{ct>,m) = 2 m g{^m). 



(14) 



Let us show that we can reach an exponential violation 
with a simple analytically manageable example. For the 
simple case of three parties (to = 3), the Mermin-Klyshko 
inequality then reads £>3 = | {-B3) | < 2, with 

B 3 =O x ®0 2 ®0' ? , + Oi®0' 2 ®O z + 0' 1 ®02®0 3 
- 0[ <g> 0' 2 ® 0' 3 (15) 

Considering a tripartite GHZ state, 

1 



GiI Z3 ) = _(|000) + |lll». 



(16) 



we have Co = c± = 2 1 / 2 and c r >2 = 0, so that 0(0, 3) = 
0i )O (0,3) = (27r)~ 3/2 cos(0) and Eq. (QJ]) becomes 

E{4>, 3) = 2 3 0(0, 3) = (2/tt) 3 / 2 cos(0). (17) 

The GHZ-like angles (6»i = 0, 6 2 = tt/6, 3 = 2tt/6, 
and 8' i — 8i + 7r/2), give the maximum violation of the 
inequality, namely 



B 3 = |3£(tt,3)- £(0,3)| =4 



3/2 



2.032. (18) 



Now, consider the multipartite generalization of the GHZ 
state 



IGHZ, 



1 



V2 



G0...0) + 11...1)), 



(19) 
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Eq. ([M]) becomes 

E(<f),m) = (2/7r) m/2 cos</>. 



(20) 



The dependence on the angle <f> of the above correlations 
is the same as the one appearing in a standard spin-like 
test for a multipartite GHZ state, namely E((f>, to) — 
cos(<p). In that case, it is known that the choice of GHZ- 
like angles k = (-l) m+1 7r(fc-l)/(2m) and 6' k = 6 k +Tr/2, 
gives the highest value of the Bell factor, namely 2( m+1 )/ 2 
(see e.g. Ref. [Hj]). Therefore, using the same angles, the 
corresponding Bell factor reads 



B m = V2 (4/7r) m/2 , 



(21) 



giving rise to an exponential violation of local realism. 

Apart from this simple analytical example, one can use 
a numerical approach to show the exponential violation 
of local realism as formula (I14|) is easily amenable to per- 
form numerical calculations for a fixed number of parties 
to. In order to find the state \^f) (coefficients c r 's) that 
maximally violates the Mermin-Klyshko inequality, one 
has to evaluate the corresponding Bell factor B m for a 
given configuration of measuring angles. The Bell fac- 
tors are expressed in general by a linear combination of 
correlation functions given each by Eq. (|14p. with the 
prescription given in Eq. ([TJ) . 

Let us search the state which maximizes the viola- 
tion of the Mermin-Klyshko inequality for the particular 
GHZ-like choice of angles. For three modes (to = 3), 
defining the (infinite dimensional) real symmetric matrix 
B 3 as 



[S 3 ] r , s = 2 3 (3 flf . 1 ,(7r,3)- flr ,.(0 ) 3)) ! 



(22) 



with the diagonal elements being set to zero, we note 
that the Bell factor can be re-expressed as S3 = C 1 B3C, 
where the elements of the vector C are given by the co- 
efficients of the input state, i.e. [C] r — c r . Consequently, 
the maximal violation of the Mermin-Klyshko inequality 
is simply given by the maximal eigenvalue of the matrix 
-B3, while the optimal input state is determined by its cor- 
responding eigenvector. In order to perform a numerical 
analysis, one has to truncate the Hilbert space dimension 
of |\&) to some arbitrary d. For example, for d = 2 the 
optimal choice turns out to be the GHZ state p^|) . giving 
a violation of B3 = 2.032. By increasing the dimension d, 
the asymptotic violation is given by B3 ~ 2.205. In Fig. [2 
we show the coefficients c r for the optimal state |\E') in 
the case d = 20, for which the Bell factor is B3 ~ 2.204. 

The same procedure can be applied for any number of 
parties. In the case of two parties (to = 2), one recovers 
the results given by Munro in Ref. [1] provided that the 
constraint c r > is taken into account, namely B2 — 
2.076. Interestingly, a higher violation can be achieved if 
we consider negative coefficients for As an example, 
we report in Fig.fJlthe coefficients c r for the optimal state 
in the case d = 30, for which the Bell factor raises up to 
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FIG. 1: Coefficients c r for the optimal state |^) in the case 
of m = 3, d — 20 (with GHZ-like measurement angles). The 
Bell factor is B 3 ~ 2.204. 
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FIG. 2: Coefficients c r for the optimal state |*) in the case 
of m = 2, d = 30. The Bell factor is B 2 ~ 2.1. 



B2 — 2.100 (recall that the Bell factor can be written 
in this case as B2 = 3£7(0, 2) — E(34>, 2), where we have 
chosen cf> = 7r/4, as in [J]). 



III. MAXIMAL VIOLATION 

Let us now see if we can find a class of states and a 
binning strategy that allows for a maximal violation of 
the Mermin-Klyshko inequality for any number of parties 
to and using quadrature measurements. We recall that 
the maximal quantum violation of these inequalities is 
given by B™ ax = 2 (m+1 ^ 2 [12]. Inspired by the results 
of Ref. M , we introduce the state 



1 



m = _l_ ilf) ®m +e i0 lg) ®m h 



(23) 



where / is a real and even function of some quadrature 
x, while g is real and odd (thus, / and g are orthogo- 
nal). Both functions are also normalized to unity. Note 
that because f(x) is real and even, it has a real and 
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even Fourier transform f(p), while g(x) has an imagi- 
nary Fourier transform ih(p), with h(p) real and odd. 

Each party t chooses to measure one of two conjugated 
quadrature via homodyne detection, either X(0) = X or 
X(tt/2) = P, and obtains a continuous variable xt or pt 
depending on the choice of measurement setting. Sup- 
pose the fc first parties measure the X quadrature, while 
the remaining m — fc measure the conjugate quadrature 
P. The joint probability that they obtain the results 
xi, ...,x k ,Pk+i, —,Pm is given by 

V{xi, ...,x k ,Pk+i, -,Pm) = \(xi\--(x k \(pk+i\---(p m \^)\ 2 

J-(f{ Xl )...f{x k )f 2 ( Pk+ l)-~f{Pm) 

+ g 2 {xi)...g 2 {x k )h 2 (p k+1 )...h 2 (p m ) 

+ 2cos(6> + (m - fc)-) x f (x 1 )g(x 1 )...f (x k )g(x k ) 

x f(Pk+i)h(p k+ i)...f(p m )h(p m )) ■ (24) 

To exploit the parity properties of / and g, we introduce 
the root binning defined in Ref . Q . This binning depends 
on the roots of the known functions / and g. If party t 
measures the X quadrature, the result will be interpreted 
as a if the measured value xt lies in the interval 

where f(xt) and g(xt) have the same sign, and "—1" if 
their signs are opposite, i.e. we consider the following 
domains 



d: 



{x e R\f(x)g(x) > 0} 
{x e R\f(x)g(x) < 0} . 



(25) 



We can similarly define the domains D+ and D~ associ- 
ated to the measurement of the quadrature P. For the 
choice of measurement settings defined above, we can 
thus calculate 2 m probability V& corresponding to the 
observation of a given collection d of binary results. For 
example, the probability V+i +i that each party ob- 
serves a result reads 



/ dxi.. 


/ dx k / dp k+1 .. 


/ dp m 




Jd+ Jd+ 





(26) 



x V(xi, ...,x k ,Pk+i, 

We are now in the position to calculate the correlation 
function E{X±, ■ ■ ■ X kl P k+ i, ■ ■ ■ P m ). Note that since / 
and g (/ and h) are even and odd respectively, f 2 and g 2 
(f 2 and h 2 ) are even functions. Hence the first two terms 
of the right hand side of (f2"4"| are even functions also, and 
their contribution to the correlation function will vanish. 
We thus obtain the remarkably simple expression 

V k W m ' k 



E{X l 



where 



X k ,P k 



fc, -ffc+i, 



V 



Pm) 



X COS 



(m-fc)| 



(27) 



W = 



dx\f{x)g(x)\ 
dp\f(p)h(p)\ 



(28) 



Interestingly, the correlation function (|27j) only depends 
on the number of sites where X and P are measured. All 
correlation functions corresponding to k measurements 
of the X quadrature, and m — k measurements of the 
P quadrature are therefore equal. We will denote them 
E(k, m — k) to emphasize this property. 

Let us illustrate the power of this compact notation 
with an example. For m = 3, the Bell factor reads 

B 3 =\E(X 1 ,X 2 ,P 3 ) + E(P 1 ,X 2 ,X 3 ) + E(X 1 ,P 2 ,X 3 ) 
-E(P 1 ,P 2 ,P 3 )\ 
=\3E(2, 1) -£(0,3)| 

7T 

2 



=\3V^Wcos( 



W 6 cos( 



(29) 



We see that the maximal violation, i.e. £?™ aa: = 4, can be 
reached with a state such that sin(#) = ±1 and V 
\\ 1 . Although such a state is quite unrealistic, one 
can define a family of physical states that approximates 
it arbitrarily well. The corresponding / and g functions 
are trains of gaussians, and V, W — > 1 as the number of 
peaks goes to infinity. We refer the reader to Ref. Q for 
their exact analytical expression. 

Let us now generalize this result for an arbitrary to. 
First note that the Bell factor @ can be written as 

B m = -\{X m B 

-(P m B' m -i)\ (30) 

with E>i = 2 Xi and B' x = 2Px- In order to benefit from 
our compact notation, we explicitely expand the expec- 
tation values of B> m -\ and B' m _ 1 in terms of correlation 
functions 



m—1 

(-B TO _i) = a k E(k, ni 1 k) 

k=Q 
m — 1 

(B' m _ 1 )=J2 a kE(m-l-k,k) (31) 



k=0 



where the a k 's are some known coefficients. When 
to — 1 = 3 for example, we have a\ = 3, a 3 = — 1, and 
<^o = a 2 = 0. As the correlation function only depends 
on the number of X and P measurements, the average 
values of the four operators of (f3"0"|) can be easily calcu- 
lated from (_B m _i) and (_B^ n _ 1 ). Suppose B m -\ has a 
term proportional to the X\...X k P k+ \ ...P m _i operator, 
which leads to the correlation function E(k,m — 1 — fc). 
The operator X m B m -\ will thus have a term propor- 
tional to X m Xi...X k P k+ i...P m _i leading to the correla- 
tion function E(k + l,m — 1 — fc), i.e. at the level of 
corelation functions we only need to replace k by fc + 1 as 
the X quadrature is measured at one additional site. A 
similar argument for the expectation values of P m _B m _i, 
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XmB^^ and P m B' m _ x leads to 



. 771 — 1 

B m =-| a k [E(k + l,m-l-k) + E(k,m-k) 



fe=0 

+E(m — k.k) — E(m 



hk + l) 



(32) 



To maximize this expression, we note that for two and 
three parties the maximum violation is reached for a state 
with V = W = 1. We thus make the reasonable hypoth- 
esis that it remains true for an arbitrary to. Recall that 
we know how to choose / and g such as to reach these 
values. When V — W — 1, we have 



which is the known maximal bound imposed by quan- 
tum mechanics. Remarkably, the state \ty) defined in 
([23| combined with homodyne detection and a binning 
strategy called root binning allows for a maximal viola- 
tion of the MK inequality. This result shows that even 
if the binning process discretizing the result of the ho- 
modyne detection discards some information, it does not 
prevent to maximally violate tests of local realism based 
on discrete variables. 



E{k, m — k) = cos 
and Eq. (f3"Tj) becomes 



(B m -i 



(B' m -x) 



m — 1 

&k cos 



k=0 
m — 1 

E 

k=0 



ak cos 



(m- fe)| 



+ (to — 1 — k)- t 



+ k 2. 



(33) 



IV. EFFECT OF NOISE 



(34) 



Introducing Equation (|33| into Eq. (|32|) and using some 
well known trigonometric formulae, the Bell factor sim- 
plifies to 



cos 

m — 1 



(o + m^j + sin (j) + to— 



x 



^ a k cos 

k=0 



(m-2k)- 



(35) 



Maximizing the violation of local realism boils down to 
finding the optimal phase 9 m such that the first factor of 
the right hand side is maximum. This term achieves its 
maximum of \/2 for a value of the phase 



We also note that 



(to - 2k)- 



= (1 — m) — 



t _i + (to - 1 — k)\ 



(36) 



(37) 



hence the maximum value of the Bell factor can be finally 
written as 



V2 



m — 1 

E 

fc=0 



ak cos 



t _i + (m - 1 - 



= V2 B, 



■max 
m—1 



(38) 
(39) 



where we have identified the summation of Eq. (|3"8|) with 
Eq. |34|) at the optimal angle 9 m -i- Introducing now the 
maximum value obtained in Ref. [5[ for the two party 
case, B r 2 lax — 2^/2, we obtain by recursion 



(40) 



In Sections [0] and IIII| we proved that the search of 
loophole-free Bell tests might benefit from an increased 
number of parties involved in the experiment. The signa- 
ture of this improvement lies in the exponential increase 
of the Bell factor with the number of parties to. How- 
ever, what makes a Bell test challenging in practice is 
not the magnitude of the violation, but rather the in- 
evitable noise associated with any real experiment. In 
many cases, these imperfections are sufficient to hide the 
non-local correlations that one tries to observe. When 
the number of parties involved in a Bell test increases, 
so does the fragility of the state used in the experiment. 
The risk is thus to rescale the violation in such a way 
that no benefit of a larger to is witnessed in practice. 
One can therefore correctly argue that an increased vio- 
lation of local realism is only significative if accompanied 
by a comparable improvement of the robustness to noise 
of the test. After all, Bell tests have to be verified in a 
lab, not on paper. 

In a discrete variable setting, the question of the toler- 
ance to noise of a Bell test is often investigated introduc- 
ing the noise fraction [l7l ]. The noise fraction quantifies 
the maximum amount of depolarizing noise one can add 
to an entangled state and still detect non-local corre- 
lations. The depolarizing noise is characterized by the 
state 1/d, where d is the dimension of the Hilbert space. 
However, in the continuous variable regime, the noise 
model underlying the noise fraction is irrelevant. Even 
if we deal with a finite number of photons, such as with 
the truncated photon number correlated states of Sec- 
tion [HJ the Hilbert space spanned by the eigenstates of 
the quadrature operators that are measured remains in- 
finite dimensional. Hence, the appropriate basis is the 
infinite photon number bases and operators proportional 
to the identity 1 have no physical meaning. To adopt an 
objective measure of the magnitude of the violation of lo- 
cal realism, we must therefore introduce a relevant noise 
model. In this Section, we will consider a probabilistic 
erasure: with probability p, the mode of a random party 
is erased; otherwise it is untouched. This noise acts inde- 
pendently on each mode and transforms an initial state 
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|*) into 

p=(l-p) m \9)(9\ 

m 

+p(i- P r- i {j2Tr t mm®\o)t(o\} m 

t=i 

+ ...+p m |0)i(0|(8)...(g) |0) m (0| (42) 

Such a probabilistic erasure is known to appear in, e.g. 
atmospheric transmissions, and has recently been studied 
in Refs. [3 

Let us first consider the photon number correlated 
states of Section [Til an d concentrate on the second term 
of Eq. (|4Tj) . Each element of the sum corresponds to the 
erasure of one of the subsystems. So, suppose for exam- 
ple that the state of subsystem m has been erased and 
replaced by vacuum while distributing The corre- 

sponding state shared between the m parties reads 

Pi, TO = Tr tn (|*)(*|)® |0) m (0| (43) 

CO 

= M 2 l r M r l ® ••■ ® l r )™-i( r l ® \°)™{o\ 

r=0 

This state is diagonal in the photon number bases, hence 
the results of all possible measurements are equiprobablc, 
i.e., V6i, ...,6 m , P(x(6i), x(6 m )) — est, and all corre- 
lation coefficients vanish. This was to be expected from 
photon correlated states as their entanglement is truly 
m-partite; tracing out one subsystem makes the state 
become separable. Thus, this property also holds for the 
other noisy terms of (j4Tj) . so that only the erasure- free 
term will contribute to the Bell factor. We obtain 

B p = {l-p) m B m (44) 

To illustrate this result, consider the m-partite GHZ 
state, Eq. (fl9]) . The noisy Bell factor reads B p = (1 — 
p) m V2(±) m/2 , hence the maximum probability of erasure 
Pmax such that non-local correlations can be detected is 

JW = 1-^2V^ (45) 

This value exponentially tends towards 1 — \pn : /2 as m 
goes to infinity, and we observe the desired increased ro- 
bustness to noise as m becomes large. 

Finally, consider now the states of Section IIIII First 
note that (f\g) = 0, hence in the {|/) ; \g}} bases, these 
states look like GHZ states. We thus expect their robust- 
ness to noise to behave like the photon correlated number 
states. Indeed, going through the calculation, one finds 
that for every noisy term of (|4ip , and for every choice 
of measurement setting the probability is an even func- 
tion of the results, e.g. if party m looses its mode, the 
probability P(xi, x m ) to obtain x\,...,x m given that 
Xi, ...,X m is measured is an even function. As a result, 
none of the noisy terms contribute to the correlation coef- 
ficients, and the noisy Bell factor is again given by (|4"4")l as 
expected. As a conclusion, with respect to probabilistic 
erasure, more parties means more robustness. 



V. A CANDIDATE 3-MODE STATE 

Let us now see how the results of Sec. IIIII can be ex- 
ploited in order to prepare multipartite states such that, 
on the one hand, they exhibit a significantly high viola- 
tion of local realism and, on the other hand, they may be 
generated with near future technology. In particular, we 
will focus on a class of three-party states whose gener- 
ation involves four optical "Schrodinger-cat" states, i.e., 
four single mode superpositions of coherent states [8( . 

Before proceeding let us make a remark. The main 
goal, as said, is now to exploit the gain expected in the 
amount of violation of locality when the number of in- 
volved parties is increased. To this aim, a natural can- 
didate may have been the generalization of the photon- 
subtracted state discussed in Refs. [f| 0- Unfortunately, 
in the case of three parties (and, actually, in the case 
of any odd number of parties) the generalization of the 
strategy adopted in Refs. d, 0] is not effective. This is 
due to symmetry reasons and the use of the sign binning. 
Specifically, on the one hand, the joint probability distri- 
bution P(x(9\) 1 x(9 m )) of any number of quadratures 
is an even function under the exchange of the arguments, 
on the other hand, the sign binning introduces an odd 
function (for odd number of parties) in the integration 
that has to be performed in order to obtain any correla- 
tion function. As a consequence, all the correlators are 
zero giving trivially no violation. 

Let us now see how the increase in the number of par- 
ties involved in a Bell test can be exploited using a dif- 
ferent approach, based on the results given in Sec. IIIII 
In Ref . [5j , the authors propose to use a superposition of 
Gaussians to implement the functions / and g of Eq. (|23[) . 
In particular, in the case of only two Gaussians they con- 
sidered the family of states defined by 

f{x)={x\[c+{\a) + \-a))\, (46a) 
g(x)={x\[c-(\a)-\-a))}, (46b) 

where c\ = 1/[2(1 ± e~ 2 l a l 2 )]. One can then calculate 
the corresponding V and W coefficients using Eq. (JS5J), 
which for large amplitudes \a\ — > oo gives V = 1 and 
W ~ 0.64. As noticed in Ref. @, in the case of two 
parties no violation is possible, i.e 2?2 — 1.90. However, 
as can be seen in Eq. (j2"5)) , one can achieve a violation of 
£>3 ~ 2.23 already for three parties. The corresponding 
state reads 

|¥ 3 ) = -L [ c 3 (| Q) _ |_ Q)) 83 + c 3 +(|a) + |_ a)) »3] _ 

' " (47) 
where we have put 8 = 0. 

Since the maximum violation is achieved for large am- 
plitudes, we can consider the following simpler state: 

|* 3 ) = c'(\a, a, a) + \a, —a, —a) + 

|— a, a, — a) + |— a, — a, a)) , (48) 




|a| 

FIG. 3: Bell factor B3 for the state |^ 3 ) as a function of the 
amplitude a [see Eq. (|48|l ]. 



S 




FIG. 4: Schematic of a possible way to conditionally generate 
the state |* 3 ) in Eq. (|48|) : (|SCS)) superposition of coher- 
ent states [see Eq. (|50[) ]; (BS) balanced beam splitter; (D) 
homodyne detector. 



preparation is experimentally demanding. However, their 
generation in traveling light modes has been recently 
reported by several groups [8j|. Thus, one may envis- 
age that, in the future, the implementation of the whole 
scheme of Fig. [4] will be possible. 



which coincides with {^3} for |ev| — > 00 [where we defined 
c' 2 = 1/[4(1 + 3e- 4 H 2 )]. In order to obtain the Bell 
factor B3 corresponding to such a state, we calculated 
the probabilities of the binned outcomes using Eq. (146)) 
to define the roots, in combination with Eqs. (|24p and 
Specifically, we considered the domains and 
Dp inherited by the state ^3): 

D+ = {x e R\x > 0} (49a) 

D~ = {x e R\x < 0} (49b) 

D+ = {p e R| - cos(pa) sin(pa) > 0} (49c) 

D- = {p e K| - cos(pa) sin(pa) < 0} . (49d) 

The Bell coefficient B3 calculated with such a proce- 
dure is shown in Fig. [3] One can see that for amplitudes 
as small as |cv| ~ 1.1, the state \^' 3 ) already gives values 
above the local bound. We note that in this regime of 
small amplitudes, \^' 3 ) 7^ ^3) and the domains defined 
in Eqs. (|49p might be non-optimal. As a is increased, 
a violation around 10% of the MK inequality is rapidly 
achieved. 

Now let us describe how the state \^>' 3 ) may be condi- 
tionally generated by using linear optics and superposi- 
tions of coherent states (SCS) of this form: 

|SCS) = c+(|a) + \-a)). (50) 

Consider the scheme depicted in Fig. |4j Two copies of 
the state |SCS) in mode ao and a\ are mixed in a bal- 
anced beam splitter. The same action is performed on 
modes 0,2 and 0,3. Then modes a[ and a' 2 , as well as 
a' and a' 3 are respectively mixed by means of two other 
beam splitters. As a last step mode o/q is measured via 
a homodyne detector. It is straightforward to show that, 
when the measurement outcome —a is obtained, then 
the conditional state of the remaining modes coincides 
approximately with |5 , 3 ). 

The main issue in implementing the scheme above 
concerns the generation of the four states |SCS), whose 



VI. CONCLUSIONS 

In this paper, we have investigated the violation of 
local realism in Bell tests based on homodyne mea- 
surements performed on multipartite continuous- variable 
states. We have proven that the Mermin-Klyshko in- 
equality supplemented with sign binning can be violated 
by an amount which grows exponentially with the num- 
ber of parties. Furthermore, we have shown that it is pos- 
sible to attain the maximal violation allowed by quantum 
mechanics by tuning the state and the binning strategy 
appropriately. The benefit of this multipartite approach 
was then shown to be effective in practice by analyzing 
the increased robustness to noise of the Bell tests with 
the number of modes. 

It is worthwhile noting that our results are not based 
on a direct mapping from continuous to discrete vari- 
ables. Such a procedure would in fact force a particular 
measurement to be used in the locality test, as the anal- 
ysis of Ref. [IH shows. We instead consider homodyne 
measurement because of its high detection efficiency. In 
this case, the possibility to obtain a maximal violation 
of the Mermin-Klyshko inequality is non-trivial, since 
(i) homodyne measurement represents a small subset of 
all possible measurements and (ii) the binning procedure 
may cause an irreversible loss of information. As it turns 
out, however, such a loss is not detrimental if suitable 
binning procedures are used, properly adapted to the 
states under investigation. 

The experimental significance of our findings lies in 
that one can effectively benefit from the increased viola- 
tion of locality in the considered multipartite continuous- 
variable scenario. In general, this allows for a greater 
freedom in the search for feasible multipartite states that 
can be used to test nonlocality in an actual experiment 
suffering from noise. As an illustration, inspired by the 
analysis of the states that maximally violate the local 
bound, we have proposed an approximation of the latter 
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state which may be experimentally realized in a reason- 
able future. This state involves three parties and gives 
a violation of around 10% of the local bound. As a per- 
spective of this work, it is expected that other states of 
light with similarly high violation and suitably adapted 
to experimental constraints might be found, inspired by 
our results. 
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